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Abstract: We study lump solutions in nonlocal toy models and their cosmological appli- 
cations. These models are motivated by a description of D-brane decay within string field 
theory framework. In order to find cosmological solutions we use the simplest local approx- 
imation keeping only second derivative terms in nonlocal dynamics. We study validity of 
this approximation in flat background where time lump solutions can be written explicitly 
and work out the validity of this approximation. Finally, we show that our models at large 
time exhibit the phantom behaviour similar to the case of the string kink. 
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1. Introduction 

In the last years strings and D-branes attracted cosmological applications related with the 
cosmological acceleration [|J-||. The observations suggest that the Universe is presently 
accelerating [||, |?]]. It is believed that a new physics is required to explain this phenomenon. 
The bulk of energy density in the Universe is gravitationally repulsive and can be repre- 
sented as an unknown form of energy (dark energy) with negative pressure and 

negative equation of state parameter w = p/p. The current experimental data show 
that w lies in the range —1.61 < w < —0.78 [j|, [l(], [|. There are various theoretical models 
- a quintessence scalar field [O], Dirac-Born-Infeld action |H]-p|, and other which are able 
to describe the case w > — 1. 

The most exciting possibility would be the case w < — 1 (see [11|-[23],|5| and references 



therein). There are several phenomenological models describing this phantom Universe 
fllfl , EH] . In these models week energy conditions p > 0,p + p > are violated and most 



of them are unstable and strange phenomena related with unstability appear [17]. There 



are also models related with modified gravity 24], but there is a problem of getting these 
models basing on fundamental principles. 
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Recently it was shown that the string field theory (SFT) description of D-brane decay 
leads at large times to an effective phantom model @. As a result of this decay the spectrum 
becomes stable and the model does not suffer from instability problems. D-brane decay in 
non-flat metric within the SFT framework is described by the string tachyon action || 

S = J ^~gd\ ^-^R+^(-^-g^d^d u <P + ^ 2 (1.1) 

Here it is supposed that we deal with 3-dimensional D-brane, G is a 4-dimensional gravita- 
tional constant, a' is a string tension, is a metric for dimensionless coordinates, R is a 
curvature, go is a dimensionless string coupling constant, fields (ft an d 3> are dimensionless 
fields related via 

(3 and k 2 are constants whose values depend on the type of a string under consideration. 
The form of the potential also depends on the type of a string. 



$ = exp(^D g )0, D g = -—d^yf^grdv, (1.2) 



Action (1.1) in the Friedmann background 

ds 2 = -dt 2 + a 2 (t){dx\ + dx\ + dxj), (1.3) 

gives Friedman equations and the equation for the field <!> = &(t) which in turn reads as 

(k 2 V + l) e - /M3 $ = U'($), (1.4) 

where V = -d} - 3H(t)d t , H = d/a is a Hubble parameter and U'(&) = dU/d&. 
To study solutions of equation (|0|) the following approximation was used 

(- K 2 + /3^)($ + 3F6) = V'($), (1.5) 
V(*) = U(*)- (1.6) 

aCq is cl parameter close to 1. This second order derivatives approximation is called me- 
chanical or local approximation. In a rather rough approximation one can take Kq ~ 1 (a 
so-called direct mechanical approximation). 

For [)kq > k 2 we get ghost kinetic term in ( |l.5|) . For this case after time rescaling we 



can take /3kq — k 2 = 1. This case corresponds to a phantom model. The dynamics of the 
phantom model in the flat background reproduces dynamics of the usual particle in the 
overturn potential |35|]. The equation of state parameter w for the phantom in a general 
potential can be represented as 

1 6 2 

t " = - 1 --K^' (L7) 

where m 2 is a dimensionless parameter related with the Plank mass M 2 , m 2 = M 2 o! = 
a'/8irG. 

For a cubic potential U the action (|1.1|) is the curved version of the action for the 



tachyon field in the Witten SFT |25] at the lowest order in the level truncation scheme 
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1 26, 27]. It is supposed that we deal with 3-dimensional D-brane in 26-dimensional space- 
time and the volume of compactified 22-dimensional space in ( |1.1| ) is omitted. 

For a quartic U the action ( |1.1[ ) is the curved version of the action for the tachyon 
field which was obtained in an approximation of slow varying auxiliary field |^] from 
cubic fermionic SFT p8[ . This action in the flat case has been obtained for the case when 
the string field contains only two fields: the tachyon in the GSO(-) sector 1 and a lowest 
auxiliary field in the GSO(+) sector. Integration over this auxiliary field produces a fourth 
degree of the tachyon self interaction. 

The existence of the kink solution for the flat version of equation ( |l.4| ) with k = has 



been shown numerically in [37] and has been proved in [ 39 1 . A validity of the approximation 
( |1.5| ) for this kink solution has been studied in [38]. 

In the case of the flat metric and k = the kink solution of (11.41) is a monotonic function 



interpolating between two different vacuum solutions. However not only monotonic kink 
solutions appear in the dynamics of string tachyons. For k / 0, k 2 < (3kq and H = 
there are non-monotonic kink solutions of (1.4) and these solutions have bounce points. In 
the open-closed string model [40] a bounce point also appears jllj]. Typical examples of 
solutions with one bounce point are lump solutions. In particular, such type of solutions 



has been recently found in [42] 



In this paper we consider two SFT inspired models with the action (1.1). For the first 
model 



U 



k + 1 



fc+i 



k < 1, 



J1.8) 



and for the second model the potential is a special polynomial 



U 



ce r . 



N 

^ n+1 

n=l 



n+1 



Distinguished feature of these models is that they have explicit lump solutions for k = in 
flat background. We study a validity of the mechanical approximation to these solutions. 
For this purpose we compare essential physical characteristics such as energy and pressure 
of the nonlocal systems with their mechanical analogues. We find that under special 
conditions these characteristics in both models are rather similar. Also in both cases 
it is possible to find small explicit modifications of the mechanical potentials to solve 
corresponding local equations in the Friedmann metric. 

The paper is organized as follows. In Section 2 we consider the equation of motion 



where < k < 1, a > 0. There is the lump solution to this equation 

*(*) = 



4f3k ' 



(1.10) 



(1.11) 



The energy of the solution (1.11) is equal to zero 



1 To produce nonBPS-branes the GSO(-) sector is involved 



(see reviews |3l| 



for details). 
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In Section 2.2 we compare ( PH) with zero energy periodic trajectories in the direct 



mechanical approximation of (|1.5[ ). The period Tfc becomes large when k — > 1. We show 
that for time t < Tfc/2 the zero-energy periodic solutions to the mechanical problem rather 
well approximate the lump solution ( jl . 1 1 ) . 



In Section 2.3 we calculate explicitly the pressure on the solution (1.11) and in Section 



2.4 compare it with the pressure of the zero-energy periodic solution of the mechanical 
problem for \t\ < T^. We find that the mechanical approximation represents a behavior 
of the pressure of the non local problem with a good degree of accuracy for k closed to 1, 
k < 1. We also compare with $ = e' 3 / 29 * $ and conclude that a comparison of $ and 
<]? gives a rather good quantitative criterium for a validity of the local approximation. A 
precise condition of a validity of the local approximation is a condition that the magnitude 
of <l - $ - §<9 2 <i> is small. 

In Section 3.1 we consider a modified integral equation with the constant friction 

e ^ t +hd t ^ h = Skg(t)<5> k h . (1.12) 

and special time depending coupling constant g(t) = e —2h(i-k)t_ (1.12) has an explicit 
lump solution 

h 



^ h (t) = e^ t+Bt , B = —(-!- Vl + 16e/? 2 ), (1-13) 

We see that the friction makes behavior of the solutions of the nonlocal equation more 
closed to that of the corresponding mechanical problem. This give us a possibility to use 
this approximation for Friedmann equations. 

In Section 3.2 we study numerically the corresponding Friedmann equations and find 
that there is a regime, i.e. suitable initial conditions, that represents an acceleration with 
w < -1. 

In Section 3.3 we reconstruct explicitly the potential for which ( 1.11| ) is an exact 
solution of the Friedman equations. This potential has the form 



where 



and 5V($>) is of order of 1/m 2 . 



V($)=V (*) + V 1 ($), (1.14) 
V ($) = -2£$ 2 log$, (1.15) 



Vi($) = ^ ty^ e r/(V-21og$) - 4<D 2 V-log^) 2 - 



P" 

V 

In Section 4 we search for lump solutions of the following equation 

A' 



S = J]a n $ n (1.16) 



n=l 



where a n are constant. It is know that a direct numerical search for lumps solutions 
is a rather difficult problem. There are also nonexistence theorems for a wide class of 
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potentials [36, 39]. It is also interesting to find lump as well as kink solutions in vacuum 
SFT fH, g]. 

Here we accept the following strategy. We start from a given function 



(fo(t) = sech 2 (i) 



(1.17) 



and find constant a n from a requirement that the function (1.17) is a solution to equation 
(1.16) with a small discrepancy for fixed N. We mean that the discrepancy has a small 
L2-norm. We consider the cases N = 3, ...,14. Qualitatively one can say that lowest 
coefficients a n show a stable picture when N — > oo. With the obtained a n in Section 4.2 
we find numerical solutions to ( |1.16| ) with zero energy using an analog of Freedman method 
(see 1 41] and references therein). 



Then in Section 4.3 we study a mechanical approximation to equation (|1.16|) . It occurs 
that all potentials V have local maximum at $ = and are equal to zero at ^oa^ &on < 1 
and ^oiv-i — &on- As has been mentioned above, the lowest coefficients a n are stabilized 
as N — > oo. For the zero-energy solutions of the mechanical problem with the initial data 
<3?(0) = one has \&(t)\ < 1 and for such $ only the low-order coefficients dominate in 
the potential. Therefore, in spite of global change of the behavior of potentials with N 
increasing, the form of potentials near the region, where the particle with the above initial 
data moves, does not change. For small time trajectories of the mechanical problems 
with above mentioned initial data do not reproduce the lump (1.17) since sech 2 (i) has 
derivatives of the same order as itself. However for large time when trajectories are closed 
to an attractor point $ = they are also closed to the lump. 

The profile of the potential ( |1.6| ) with found a n in the region < < 1 is rather 
closed to the profile of the potential 



y ($) = 2 (1 - $) $ 2 



(1.18) 



for which the function ( |1.17 ) solves the mechanical problem with the lump boundary condi- 
tions. In Section 4.5 we show that engaging of the Friedmann metric modifies the potential 



(1.18) and one has to add 



Vi=- 



75m| 



(-l + $) 3 (2 + 3$) s 



to ( 1.18| ) in order the function ( 1.17j ) to be a solution of the Friedmann equations. 



2. Gaussian Lump 



2.1 Action and equation of motion 

In this section we consider the action 



S 



k + 1 



fc+i 



(2.1) 
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with an unusual non-polynomial interaction 

< k < 1 

and where □ = — d 2 + V 2 , 

$(x) = e2 n (f)(x), a > 0, 
The equation of motion for this action is 

For space homogeneous configurations $ = it takes the form 

here and below d = d^. 

Equation (2.4) has the following solution 

b_ t 2 1 - k 

$(t) = e 4 /3 , 6 = — — 
k 

The relation between <& and </> is 

$(t) = e -§ a2 0(i) 

and one can get from $ using the smoothing integral representation 



i r it-*') 2 



e ^ <f>(t')dt' 



Substituting (|2.6|) into (|2.8|) one gets smoothed solution 



0(t) 



\/2 + & 



g 2/3(2+6)' 




(2.2) 
(2.3) 

(2.4) 
(2.5) 

(2.6) 
(2.7) 

(2.8) 
(2.9) 




\\ -6 -4 -2 2 4 6 2) 

Fig. 1: Solutions of integral equation (|]|): 1) for = \, k = 0.2; 2) for = \, k = 0.8. 
The firm line is and the dashing line is <f>. 



From ([O]) and ( |2,9| ) we see that and cf) are related as 

2k 

<j)(t) = y/c$(y/ct), where c 



l + k 



(2.10) 



When k — ► 1, c — ► 1 and therefore = $ 2 and when A; — > 0, c — > and = 0. This 
is illustrated on Fig. [I] where solutions of ( |2.5| ) and corresponding smoothed fields 4> are 
presented for different values of k. 

2 Let us note that in the more complicated model E] (ca = ^r) it is known, that physical (analog of <f>) 
and tilded (analog of $) fields coincide with very high precision 
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2.2 Mechanical approximation 

Let us consider the mechanical approximation for equation ( |2.5| ), 

m) = -vm, (2.11) 

the parameter (3 plays a role of mass of the particle. The potential V t is an overturn 



version of the potential (1.6) for U given by (|TT 
V ot (<p)=-V( V ), V(<p) 



k + 1 



fe+i 



1 2 

2* 



(2.12) 



To have a real potential V in ( |2.12| ) let us assume that k = 1 — 26 and present the 
corresponding square in U as (<3? 2 )1/(1 — 5). In this case the potential V^($) is an even 
function and this function is equal to zero at $ = and has beak with singular second 
derivative $ = 0, see Fig.^.a. The particle with the zero energy oscillates between the point 
(p = and the bounce at <p ma x = (f^pf) 1 " 7 "- The zero-energy trajectories are presented in 
the fig.^.b and the exact lumps ( |2.6D are presented on fig.||.c for different k. 




c) 

Fig. 2: a) The potential V for the k= 0.2 (dashing line) and for k=0.8 (firm line), b) The 
numerical solution of equation ( p. 11 ) for k=0.2 (dashing line) and for k = 0.8 (firm line), 
c) Exact solution to integral equation ( |2.5[ ) for = 1/4 and k= 0.2 (dashing line) and 
k=0.8 (firm line). 



The periods of these motions are 



dip 



/o V—ZVot 

In the table below we enumerate values of T for different k and (3 = 1/4. 



(2.13) 



k 


1/5 


1/3 


1/2 


2/3 


8/9 


9/10 


T 


5tt/4 


3?r/2 


2?r 


3vr 


9tt 


10vr 



The solutions of the integral equations (2.5) (see Fig. ||.c) and their mechanical analogs 
(2.11) describe different physical behavior. Namely, the solution fl2.6| ) of the integral equa- 
tion starts from zero go above the bounce point (p ma x then rises $ = 1 and returns to 
zero. Thus the solution of the integral equation does not bounded by the bounce point, 
whereas the zero-energy solution of the mechanical model presents a periodic motion be- 
tween the bounce point and zero. The period depends on k and as k goes to 1 the period 
becomes infinite, i.e. k = 1 corresponds to the free motion. Therefore, for small k « 1 
the mechanical approximation does not work, while it works for t < for k goes to 1. 
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2.3 Energy and pressure 

In this subsection we consider the energy and the pressure of the nonlocal problem 
and compare them with energy and pressure of the corresponding mechanical problem 

(2.14) 



( 2.11 ). Equation ( |2.4j ) has the conserved energy (compare with [35]) 

E = E p + E nil + E n i 2 , 



where 



E n u 



E„ 



2^ jfe + 1 



e 2' 



(2.15) 
(2.16) 



and 



2 Jo 



(2.17) 



One can see that this energy conserves. Indeed 
dE(t) 



dt 



+ \fk$ k d$ + - [\e~2 pd2 Vk$ k )l)(e2 pd2 d<S>)dp 
2 Jo 



+ y/k$ k d& - Vk<f> k e~2 d2 dd> = 
The pressure in the model ( |2.1[) has the form 

P = E n i 2 



+ e 2 



0. 



E p — E nil 



(2.18) 



Substituting in this formula expressions for E n \ x and E n i 2 given by ( gig ) and ( gig ) one 

gets 

P = -p (e( 2 - p ^ d2 d<S>){e2 pa2 d§)dp (2.19) 
Jo 

Substituting in ( |2.19|) the explicit form of the solution (|2.6D we get 



P 



2bH 



.2+2 rl 



dp 



exp(- 



6(2 + b) 



(2.20) 



A) [(2 -b(p-2))(2 + ~"" V /3 4 + 46-6 2 (p-2)//' 

From formula ( 2.20| ) we see that the pressure is negative and for large time goes to 



zero. 



2.4 Compression of pressure for nonlocal and local problems 

Comparing expressions ( |2. 14 ) and ( |2.1§| ) with standard expressions of the energy and pres- 
sure for a scalar field, E v = ^ip 2 + V((p), P^ = \p 2 — V(tp), one can say that E n \ x and E n i 2 
play roles of extra terms in the nonlocal "kinetic' energy £k and 'potential' energy Ep, i.e. 



E — £k + £ P 
P — Sfc £p 



(2.21) 
(2.22) 
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and 



E, 



III 2 



£ P — E p + E n i x 



(2.23) 
(2.24) 



where E p and E n i x , E n \ 2 are given by ( [2.15 - 2.17 ). 

On the solution ( |2.6j ) we can write explicitly E p , E n \ x and E n i 2 as functions of time. 
These functions are presented on Fig|3| for particular value of parameters (3 = 1/4, k = 0.8 
and k = 0.2. Time dependence of £ p and £k is also presented on Fig. || 
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Fig. 3: a) The potential energy E p ; b) The nonlocal terms E n \ x (firm line) and E n i 2 (dashing 
line); c) £ p (dashing line) and £k (firm line). la,lb and lc correspond to k = 0.8, f3 = 1/4, 
and 2a, 2b and 2c to k = 0.2, 13 = 1/4 



From Fig.^.c we see that £ p > and £k < for < t < oo. Since the energy is 
conserved and equal to and the pressure is < 0. 
Representing (|2.16| ) and (|2.17|) as 

Enh = I [\e^-^ a2 ^){e^ 92 d 2 ^)dp (2.25) 
2 Jo 

and 

Em 2 = -~ / {{e^ l ~^ 92 d^)){e^ d2 d^)dp. (2.26) 
1 Jo 

we can see that in the approximation neglecting high order derivatives one has 

Em 2 ~ -^ 2 , (2.27) 

E nh ~ ^d>$, (2.28) 

The discrepancy in these approximations are illustrated on the fig.|j. We see that the 
discrepancy becomes smaller when k — > 1. 



Fig. 4: 

E n [ x and the 



On the left pictures (la and 2a) the dashing lines present the exact expressions 
sent the approximations to E n \ x . On the right pictures (lb and 
ent the exact expressions E n i 2 and the firm lines present the 
On la and lb k = 0.8, = 1/4 and on 2a and 2b k = 0.2, 



2b) 

approximations to E n \ 2 
P = 1/4- 



In this approximation 



2 



(2.29) 



p 2 r jfc + l 2 



Taking into account that in this approximation 



P. 



_i$ 2 _ — cj>$ 
2 2 



we get 



f 



2 fc + 1 



(2.30) 



(2.31) 



(2.32) 



So we see that the approximated kinetic term ( 2.2S| ) and potential £ p ( p. 32 ) coincide up 
to the sign with the kinetic and potential energy for the mechanical problem ( 2.11| ). 

On Fig|| the kinetic and potential parts of the energy for the local problem ( 2.1l| ) 



as functions of time are compared with and £ n calculated on the solution (2J3). We 
see that for k — > 1 the kinetic and potential parts for local and nonlocal problems almost 
coincide for time < t < T^. 
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0.015- 
0.01 
. 005 



-0.005 
-0.01 
-0.015 




b) 



Fig. 5: Local kinetic (dashing positive line) and potential (dashing negative line) parts of 
the energy on zero-energy solutions of eq.( |2.11| ) as functions of time. Nonlocal kinetic (firm 
negative line) and potential parts (firm positive line) of the energy on the solutions (|2.6|); 
the picture for k=0.8 (a) and k=0.2 (b). 



3. Perturbations of Gaussian Lump 



3.1 Perturbation by the friction 



Let us consider the integral equation with the constant friction and time depending coupling 
constant 



(3.1) 



For the case of a special time-depending coupling constant, 



9® 



-2h(l-k)t 



there is the following lump solution 



* ft (t) = e ^ 



±t 2 +Bt 



B = —(-l-^TTWb) 



(3.2) 



Using the following integral representation 



m+hd t 



*(«) 



40 



-$(t')dt'. 



(3.3) 



one can check that (^) solves ([O]) 

Let us consider the mechanical problem with a friction that approximates equation 

0<p(t) + fup(t) + <p(t) = g(t)Vkcp(t)\ (3.4) 

The particle with the zero energy starts from the point ip m ax = (i^x) 1 ~ fe and moves to 
ipo = 0, but do not reach it. The corresponding trajectories are presented in the fig.§.b for 
different k and on fig.||.c are presented the exact lumps (2.6) for different k. 
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Fig. 6: a) The numerical solution of equation Q3.4| ) for k=0.2 (down solid line) and for 
k = 0.8 (upper solid line) and exact solution to integral equation fl2,5|) for f3 = 1/4 and 
k= 0.2 (down dot line) and k=0.8 (upper dot line), h=0.5; b) the numerical solution of 
equation (|3.4[) for k=0.2 (upper solid line) and for k = 0.8 (down solid line) and exact 
solution to integral equation fl2.5| ) for (3 = 1/4 and k= 0.2 (down dot line) and k=0.8 
(upper dot line), h=0.1 

We see that the friction makes behavior of the solutions of the nonlocal equation more 
closed to that of the corresponding mechanical problem. 

3.2 Perturbation by non flat metric 

Let us now consider an interaction of the theory (|2.1[) with the gravity 



J s/^gd'x 



where 



$ = exp(|n fl )0, 



□ , 



fc + 1 



(3.5) 



(3.6) 



On space homogeneous configurations in the Friedmann metric ( |1.3j ) the field equation 
takes the form 

e^*(t) = Vk<S> k {t) (3.7) 
where V> = —d\ — 3H(t)dt and H(t) = a/a, a = dta. The Einstein equations have the form 



3H 2 = At p 



H 2 + 2d a = = p 



with the energy and pressure densities are given by 

P 



1, 
2 



k + 1 



where 



p = 2£ 2 

f 1 dT{e^ TV Vk<S> k )V{e-^ TV <$>) 
Jo 



P f 1 
2 



(3.8) 
(3.9) 

(3.10) 
(3.11) 

(3.12) 
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Z 2 = -t f 1 dT(e% TV dVk$ k ){e-% rV d$) 
2 Jo 



Motivated by the flat case we make the following approximation 
exp(df + 3H(t)dt)$ n(l + d? + 3H(t)d t ) $ 



The corresponding Friedmann equations have the form (3.10), ( 3.11| ) with 

p = -\v 2 + v{v), 

P = -^jp 1 ~ V{tp) 



and the equation for ip reads 
Here V is given by (2.12) 



Cp + 3Hip = V' 



(3.13) 
(3.14) 

(3.15) 
(3.16) 

(3.17) 






Fig. 7: a) $ as function on t with the initial condition </?(0) = 0.642 and </?(0) = 0; b) 
function H = H(t) as function on t; c) w = w(t) as function on t. 



An numerically solution of the system of equations ( 3.15j ) and fl3.17|) for the potential 
( 2.12 ) is presented on Fig]?} 



3.3 Exact Gaussian Lump Solution in Friedman Metric 

In section [2] we have analyzed the solution for the phantom field with fractional interaction 
in the Friedman metric. Here we shall construct the potential which leads to exact gaussian 
lump-type solution (2.6). It is known that if one has an explicit solution for the space- 
homogeneous scalar field then one can reconstruct the form of potential^, [23fl . To realize 
this one assumes that the Hubble parameter H(t) has the special form 



H(t) = W{<P{t)) (3.18) 

where W is a function called superpotential. Substituting the anzats ( |3.18| ) in the Fried- 
mann equation H = (ft 2 jlvnt we get 

2m 2 p W' = <p (3.19) 
where prime denotes derivatives in <j). For the gaussian lump solutions we have 

<p = -2<£v/-£ M. ( 3 - 20 ) 



c) 
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that gives the equation 



2M 2 W' = -2<j)y/-£ huf>, 
from which we reconstruct the superpotential 



(3.21) 



1 



(3.22) 



The potential is obtained from the second Friedmann equation, V((p) = l/2t/> 2 + 3m 2 VF 2 , 
and 



V{4>) = -2^ 2 log0 + -pL (v / 2Ter/(V-21og0) - 4</> 2 ^-log^))' 

b4mf V / 



64 m 2 



(3.23) 



0.2 0.4 0.6 0. 





-80 
-100 



a) n 1 5 5 3 w b) ~ luut ' c) 
Fig. 8: a) The overturn potential — V, b) H = H(t) as function on t; c) w = w(t) as 
function on t and the line w = — 1 



The Hubble parameter has the form 

1 



H(t) 



At£e- 2 & - v / 2<er/( v / 2^t) 



(3.24) 



The asymptotic behavior of the Hubble parameter and the scale factor are 



H(t) ~ H 



8m 2 



a(i) ~ e Hot for large t 



The state parameter has the following explicit form 



w(t) 



256^ 2 e 2 ^m 2 



3 (-At + e 2 * 2 « V2^erf(V2t V?))' 



(3.25) 



(3.26) 



The shapes of the potential, the Hubble parameter and the state parameter are presented 
at (Fig.||). It is interesting to note that the next to leading term in the asymptotic of ( 3.26j ) 
at large t 

128m 2 t 2 £e- 2 ^ 2 
1 p 



w(t) 



3tt 



(3.27) 



is proportional to m 2 meanwhile from the general formula (1.7) one can expect that the 



deviation from the limiting value -1 could be proportional to 1/m 2 . This is caused by a 
presence of a dependence of m 2 in Hq ~ and independence of 4> on m 2 . 
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4. Hyperbolic Lump 

4.1 Numerical construction of the potential 

Here we will be interested in the scalar SFT-type nonlocal models described by the following 

'1 



action 



S = d a x 



[/($) 



(4.1) 



where U is a polynomial and <3? the smoothed field given by (|2.8|). Action (4.1) leads to 
the equation of motion which on the space homogeneous configurations reads 



We assume that U has a form 

N 

tf(<l>) = ^(cl>) = V^<J> n+1 , 



n=l 



We choose a function 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



is small in L2-norm for fixed N, where N is a number of terms in the series (|4.3| ). For 
numerical calculations we use a representation of the left hand side of ( [4.2D in term of the 
integral operator 



$o(t) = sech 2 (t) 

and find the set of coefficients {a n } from a requirement that that a discrepancy 

5 N (a) = e^ *-[/'($) 



(t-T) 2 



e ^ $(t)<Zt, 



(4.6) 



and take (3 = \- We denote the R.H.S. of ( |4.2j ) as a force F 



A/ 



(4.7) 



Note that numerical calculations show that we cannot make the discrepancy equal to 
zero for final iV, this illustrate the fig.^a). 



0.01 
0.005 



2 4 6 8 10 12 14 



a) 




-4 -2 2 4 6 



b) 




Fig. 9: a) Minimum of the absolute value of the function of discrepancy 5j^(a) ( |4.5[) (in 
the logarithmic scale at pintle) depending on number of terms N is presented; b)$o(£) = 
sech 2 (i) is presented by dash line, function if[$o](^) is presented by the firm line and 
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function F a5 [ipo\ is presented by dotdashing line N = 5, coefficients «5 are collected on the 
Table; c) function S 5 (a 5 ) = K[<p ] - F a5 [<p ]. 

Note that already for N = 5 the discrepancy is rather small, this illustrate the fig|]b). 
The coefficients a for different N are presented in the Table a 



N 


ai 


Oil 


a 3 


ai 


«5 


a 6 


OL-J 


3 


1.896 


-2.2206 


1.0649 










5 


2.319 


-5.7831 


10.1891 


-9.1559 


3.1602 






6 


2.4309 


-7.5998 


19.2683 


-28.4428 


21.4779 


-6.4098 




12 


2.4628 


-8.2965 


23.6528 


-38.931 


28.6269 


0.287057 


-6.3996 


13 


2.4643 


-8.3152 


23.7106 


-38.965 


28.5816 


0.278636 


-6.3756 


14 


2.4662 


-8.3356 


23.7655 


-38.985 


28.5391 


0.256198 


-6.3702 


15 


2.4602 


-8.2969 


23.7264 


-39.024 


28.5326 


0.282346 


-6.3252 



It is interesting to note that as follows from the fig. |9[ there is a tendency of decreasing 
of the minimum value of the discrepancy when the number of terms in the series (|4.3| ) is 
increase, moreover the values of the low-order coefficients are stabilized when N is increase. 

4.2 Numerical solution of the integral equation 

We are going to solve the equation 

K[y](t)=F aN [y](t), (4-8) 

where the right hand side is given by F aN , that minimized the discrepancy for the function 
( |4.4| ) A solution of this equation we find using the following iteration procedure 

y n +i = yn + X(K[y]-F aN [y}), A = 0.01, (4.9) 

and as the zero approximation we take fl4.4| ), i.e. yo = sech 2 (t) 





400 500 



c) 



Fig. 10: a) Zero iteration yo = sech 2 (t) is depicters by firm line and five hundredth 500 by 
dashing line(iV = 5); b)The discrepancy for initial iteration (the line with amplitude 0.003) 
and for solution y(t) of equation Q4.8| )(the line with amplitude 0.001); c)The dependence 
for descrepansy ( [4.10[ ) as funcion on number of iteration when the number of terms in the 
series N = 5. 



Numerical calculations show that this iteration procedure converges. A character of 
a convergent of the iteration procedure is natural to specify by the a maximum of the 
discrepancy 

A[y n ) = max\K[y n ](t) - F aN [y n ]{t)\ (4.10) 



- 16 - 



As we see that A[y n ] becomes smaller for a large number of iteration n (when the number 
of terms in the series N = 5), however as we can see from Fig. |l^ c) at n ~ 500 A[y n ] 
goes to a "plateau" and after decreases slowly. The same situation takes place when we 
consider the case when the number of terms in the series N = 7, 14, 21, ..etc. 



4.3 Mechanical problem 

The mechanical problem that corresponds to (|4.8| ) has the form 

P<pit)=F a [<p]-<p. 
The potential for this problem is V t{<-p) 

1 



Vot(<p) = ^P 2 - U s (ip) 



(4.11) 



(4.12) 



V t(<p) is the overturn version of the potential V given by (1.6) with U given by (4.3). 






aj b) 
Fig. 11: The overturn potential (4.12) when N is equal: a) N = 5 and b) N = 14; c) the 
overturn potential — Vq 



The overturn potentials for N=5 and N=14 are presented on fig. 11. It is interesting to 



note that global forms of corresponding potentials are different for different ./V see fig. 11 



but in all considered cases there are local maximum of the potential in the zero field 
configuration and points of local minimum. 

Starting from the zero-energy bounce point, <fi(0) = (fo, where V{ipo) = 0, with zero 
velocity, <p(0) = 0, the particle moves to the local maximum of the potential at the point 
99 = during infinite time. For the initial data (f(0) = <pi < <po and (p(0) = the particle 
performs a periodic motion. For f(0) = ipo and (p(0) 7^ the motion is unbounded. 

On Fig.O,c for a comparison we present the potential — Vq. As it is known (see also 



the subsection (4.5)) the lump (4^) is a zero-energy solution of the mechanical problem 
with the phantom kinetic term and potential Vq given by ( 1.18| ). 



4.4 Energy and pressure 

In this subsection similar to the subsection 2.3 we consider the energy and the pressure of 



the nonlocal problem (4.8) and compare them with energy and pressure of the correspond- 
ing mechanical problem ( 4.1 1| ) . Equation ( |4.8| ) has the conserved energy functional 



E — E p + E n i x + E n i 2 , 



where 



E n 



(4.13) 



(4.14) 
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and 



The pressure 
has the representation ( ]4. 17| ) and 



E nh = | [\e-i pd2 U'm(e^ 92 d 2 <5>)dp 
1 Jo 

>nh = ~ (\{e-^ d2 dU'{<S>){e^ d2 d<S>)dp. 
1 Jo 

P = E n i 2 — E p — E n i x 



P = -al {e {2 - p ^ d2 d<5>){e2 pd2 d<5>)dp 




(4.15) 
(4.16) 
(4.17) 

(4.18) 



Fig. 12: a)The lump (|4.4| ) and the corresponding smothered function (lowest positive line 
in the picture) and the pressure fl4.18|) for N = 14 (negative line). 



4.5 Exact Lump Solution in Friedmann Metric 

In the full analogy with what has been done in sect. 3. 3 we can easily find a potential for 
which the function (4.4) is a solution to the Friedmann equations. It has the form 



V(cf>) =2(1-<P)<P 



2 4(-l + ^>) 3 (2 + 3 l 
75 m p 2 



(4.19) 



The shape of the potential is presented at fig, 13 .a. 




3 4 5 



Fig. 13: a)The overturn potential — V, b) H = H(t) as function on t; c) w 
function on t 

The Hubble parameter has the form 

2 (4 + cosh(2t)) sech(t) 2 tanh(t) 3 



w(t) as 



H(t) 



15 m p 2 



(4.20) 
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and its shape is presented on (Fig.|l3|,b). H(t) and the scale factor a(t) have the following 
asymptotic behavior at large t 



H(t) ~ — a(t) ~ e 1 ^ /or t oo (4.21) 

The state parameter is given by the following explicit formula 

75 ml coth(t) 4 

w ( t ) = -1 ? (4.22) 

(4 + cosh(2t)) 2 V ' 

and we see that for large times w(t) goes to —1. 
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